Classical nucleation theory in ordering alloys precipitating with Ll 2 structure. 
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By means of low-temperature expansions (LTEs), the nucleation free energy and the precipi- 
tate interface free energy are expressed as functions of the solubility limit for alloys which lead to 
the precipitation of a stoichiometric LI2 compound such as Al-Sc or Al-Zr alloys. Classical nu- 
cleation theory is then used to obtain a simple expression of the nucleation rate whose validity is 
demonstrated by a comparison with atomic simulations. LTEs also explain why simple mean-field 
approximation like the Bragg- Williams approximation fails to predict correct nucleation rates in 
such an ordering alloy. 



Since its initial formulation in 1927 by Volmer, We- 
ber and Farkas and its modification in 1935 by Becker 
and Doring the classical nucleation theory (CNT) 1,2,3 
has been a suitable tool to model the nucleation stage 
in phase transformations. The success of this theory re- 
lies on its simplicity and on the few parameters required 
to predict the nucleation rate. Recently, the use of com- 
puter simulations have allowed to assess the applicability 
of the theory for solid-solid phase transformations^^. 
Thanks to a precise control of simulation conditions, it is 
possible to get accurate estimations of CNT parameters 
and thus to make a direct comparison between theory 
predictions and quantities observed during simulations. 
One thus gains a deeper understanding of the validity of 
the different assumptions used by the CNT. 

Previous studies have shown that the capillary approx- 
imation, which CNT relies on, gives a precise description 
of cluster thermodynamics. Within this approximation, 
the free energy of a nucleus is written as the sum of a 
volume contribution, the nucleation free energy, and a 
surface contribution corresponding to the energy cost to 
create an interface between the nucleus and the solvent. 
For CNT to agree with atomic simulations, care has to 
be taken in the way these two energetic contributions 
are obtained. In particular, we have shown that one has 
to take into account short range order when calculating 
the nucleation free energy in an ordering alloy 4 -. Usual 
thermodynamic approximations, like the ideal solid so- 
lution or the Bragg- Williams approximation, cannot de- 
scribe short range order and thus can predict values of 
the cluster size distribution and of the nucleation rate 
wrong by several orders of magnitude. This is to contrast 
with more sophisticated mean-field approximations like 
the cluster variation method (CVM) which provides good 
predictions of the nucleation rate 4 *. However, an easy use 
of CNT and a clear determination of the missing ingre- 
dients in simple mean-field approximations requires an 
analytical approach which CVM cannot provide. Such 
an approach has to lead to accurate expressions of the 
CNT input parameters so as to make the theory predic- 
tive without any fitting of its parameters. 

In this Letter, we use low-temperature expansions 
(LTE) 8 ">2 to derive an analytical formulation of the nucle- 
ation free energy and the interface free energy in a binary 
system like Al-Sc or Al-Zr, i.e. a supersaturated Al-X 



solid solution leading to the nucleation of a stoichiomet- 
ric AI3X compound with the Ll 2 structure. This struc- 
ture corresponds to an ordering of the fee lattice with so- 
lute X atoms lying on one of the four cubic sublattices 8 . 
LTE are well suited to describe short range order in 
dilute solid solution and nearly stoichiometric ordered 
compound a 10 ! 11 ' 12 ! 13 ' 14 like AI3X compound. The use of 
this method in CNT framework allows to obtain a fully 
analytical modelling whose only material parameters are 
the solubility limit and the solute diffusion coefficient. 

To do so, we start from the same atomic diffusion 
model previously developed for Al-Sc-Zr system 4 ^ 5 . This 
model relies on a rigid lattice with interactions between 
first- and second nearest neighbors and uses a thermally 
activated atom- vacancy exchange mechanism to describe 
diffusion. Despite its simplicity, it has been shown to 
lead to predictions in good agreement with experimental 
dat a 15 i 16 ' 17 . Within this atomic model, atoms are con- 
strained to lie on a fee lattice and the configurations of 
a binary Al-X alloy is fully described by the solute atom 
occupation number p n with p n = 1 if the site n is occu- 
pied by a solute atom and p n = otherwise. The energy 
per atom of a given configuration of the Al-X alloys is 
then given by 



E = U M + ^- (Ux - U M ) Y\ Pn 

+ W U} AIX^2 Q--Pn)Pn 



+ , >< 2 > 



V (l-p n )Pm (1) 
* — *n,m 



where the first and second sums, respectively, run on all 
first and second nearest-neighbor pairs of sites, N s is the 
number of lattice sites, Ua\ (respectively Ux) is the en- 
ergy per atom when only Al (respectively X) atoms lie 

on the fee lattice and w^\' x and w^ix are tne nrst ana ^ 
second nearest neighbor order energies. Al-Sc and Al- 
Zr thermodynamics are characterized by the order ten- 
dency between first nearest neighbors (wjy X < 0) and 
the demixing tendency between second nearest neighbors 
(wjyx > 0). Eq. [T]is a rewriting for binary alloys of the 
atomic model developed in Refs^i^ when one neglects 
vacancy contributions. 
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The nucleation free energy entering CNT is defined by 



AG nuc (a: x ) 



- [/iAl(Xx 9 ) 



MAi(a^x)] 
i[M^)-Mx(4)]. (2) 



where /j,ai(%x) & n d Mx(^x) are the Al and X component 
chemical potentials in a solid solution of concentration 
xx, and ir x and x x the concentrations of the equilibrium 
and supersaturated solid solution. 

LTE are more easy to handle in semi-grand-canonical 
ensemble where all quantities are written as functions of 
the effective potential /i = (fi ai — Mx) /2. Definition of 
the nucleation free energy then becomes 



AG nuc (fi) = A(^) - A(fi) + - ( M eq - n) , (3) 

where /i cq is the effective potential corresponding to equi- 
librium between the Al solid solution and the AI3X LI2 
compound. We have defined in Eq. [3] the solid solution 
semi-grand-canonical free energy A = (/iai + A*x) /2 = 
F(x) + (1 — 2x)^, F(x) being the usual canonical free 
energy. 

A LTE consists in developing the partition function 
of the system around a reference state, keeping in the 
series only the excited states of lowest energies. Use of 
the linked cluster theorem^ allows then to express the 
corresponding semi-grand canonical free energy as 

A(ji) = A°{v) ~ kTj29i, n exp{-AE iin (iJ,)/kT), (4) 



where the energy of the ground state is -4°(/i) = Uai + A 4 
for the Al solid solution and A (/i) = 3/4 U Ai + 1 /4 U x + 
3 w aix + /V2 for the AI3X LI2 compound. In the sum ap- 
pearing in Eq. [4j the excited states have been gathered 
according to their energy state i and the number n of 
lattice sites involved. LTE parameters corresponding to 
the excited states with the lowest energies are given in 
Tab. U All excitation energies only involve a set of iso- 
lated atoms or in second nearest neighbor position since 
flipping two atoms at nearest neighbor position produces 
an excited state with a much higher energy. 

The solute concentration in a given phase is obtained 
by considering the derivative of the corresponding semi- 
grand-canonical free energy. For the solid solution, one 
gets 



2x(m) = \ 



1 



dA(y) \ 



^n^exp (-AE iin (p)/kT). 



(5) 



The solid solution and the Ll 2 compound are in equi- 
librium when both phases have the same semi-grand 
canonical free energy. Considering third order LTE 



TABLE I: Coefficients entering in the low temperature ex- 
pansion (Eq.Q. The first seven excited states are considered 
for the solid solution and the first three excited states for the 
AI3X LI2 compound. The effective potential is written as 
H = (U x - U A i)/2 + 6a4Yx + Sfi. 
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LI2 compound 
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(i = 3), this happens for the effective potential 



M cq = (Ux - U M )/2 + 6wj& 



-kT 



exp (SufykT) + 3 exp (-lOc^/x/ kT ) 



(2) 



— expi-Uu^/kT) 



(0) 



corresponding to the solubility 



= exp (-6u$ x /kTj + 6 exp f-lOw^/fcT 



16 exp 



-Uj^/kT 



(7) 



As these expressions have to be consistent with the ex- 
pansion of A, terms with larger exponential arguments 

(2) 

than — 12a; jy X are discarded. For equilibrium phases, one 
does not need to go further in the expansion than the 
third order. Indeed, thermodynamic properties are al- 
ready well converged as the solid solution and the LI2 
compound in equilibrium only slightly deviate from their 
respective ground states. On the other hand, Fig. Q] 
show that an expansion beyond the third order of the 
semi-grand canonical free energy A(/i) of the supersat- 
urated solid solution and of the corresponding concen- 
tration xx(fJ-) improves the convergence of the nucle- 
ation free energy. When only the first excited state 
is included in the expansion, LTE leads to the same 
value of the nucleation free energy as the ideal solid 
solution model. As more excited state are included in 
the expansion, the value deduced from LTE converges 
to the one calculated with CVM in the tetrahedron- 
octahedron approximation^. This is to contrast to the 
Bragg- Williams approximation which leads to a worse 
prediction of the nucleation free energy than the ideal 
solid solution model. 
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FIG. 1: Variation with the nominal concentration x^ T of 
the nucleation free energy AG nuc at T = 723 K obtained 
with different thermodynamic approximations: CVM, ideal 
solid solution, Bragg- Williams and low-temperature expan- 
sions (LTE) to different orders. 



So as to understand why Bragg- Williams approxima- 
tion does so bad, it is worth going back to the canonical 
ensemble. When considering only the third order LTE, 
thermodynamic quantities can be expressed as functions 
of the solid solution nominal concentration 4- In par- 
ticular, the nucleation free energy is given by 

AG^ c E (4) = fcT [q(4)-q(4 q )] 

+ sfcTexp [ q (4) 2 - q (4 q ) 2 

-^r{ln[q(4)]-ln[q(4 q )]}, (8) 

where we have defined the function 

2x 



q(x) 



1 



1 + 4x 



6 exp (2J^ x /kT 



19 



(9) 



This expression developed to first order in the concentra- 
tions 4 and ix leads to 



Comparing Eq. [10] with Eq. [TT], we see that these two 
thermodynamic approximations deviate from the ideal 
solid solution model by a distinct linear term. In the 
LTE (Eq. [TO]) . the nucleation free energy is only de- 
pending on the second nearest neighbor interaction and 
the coefficient in front of the concentration difference is 
positive. On the other hand, the Bragg- Williams ap- 
proximation (Eq. Hip incorporates both first and sec- 
ond nearest neighbor interactions into a global parameter 
lu aix = Gw^ix + 3w? lx . This leads to a linear correction 
with a coefficient which can be negative due to the os- 
cillating nature of the interactions. In particular, this is 
the case for both binary Al-Zr and Al-Sc alloys^. Bragg- 
Williams approximation thus leads to a wrong correction 
of the ideal model because it does not consider properly 
short range order. In the case of a LI2 ordered com- 
pound precipitating from a solid solution lying on a fee 
lattice, one cannot use such an approximation to calcu- 
late the nucleation free energy. On the other hand, Eq.[5] 
is a good approximation and can be used to calculate 
the nucleation free energy even when the second nearest 

,(2) 



neighbor interaction u\rv * s n °l known. Indeed, this pa- 
rameter can be deduced from the solubility limit x^ 3 by 
inverting Eq. [7J leading to the relation 



4 2 / x = -ifcTln(4«) 



kT 



cq2/3 



u cq 

3 Xx 



(12) 



This relation combined with Eq. [8] provides a powerful 

way for calculating the nucleation free energy from the 

solid solubility. 

LTE can be used too to calculate the plane interface 

free energy cioo corresponding to a [100] direction. Due 

to the inhomogeneity perpendicular to the interface, the 

main contribution arises from broken bonds and excited 

states, whose energies are lower than in bulk phases, only 

bring a small correction. One thus does not need to go 

further than the second order in the expansion. At K, 

the isotropic interface free energy a is obtained by mul- 

1/3 

tiplying crioo with the geometric factor (6 /it) corre- 
sponding to a perfect [100] facetting of the precipitates. 
For low temperatures, this is a good approximation to 
assume that the same linear relation holds between both 
quantities 4 . The isotropic interface free energy given by 
LTE is then 



1 / x v 

-fcTln 



+ hr (l + 6e 2 "A 'x/ feT ) (4 - a£) . (10) 

Doing the same development for the nucleation 
free energy calculated within the Bragg- Williams 
approximation 4 , we obtain 



(r° ) 



-fcTln f- 
4 I 1 



>m(- 

4 \ x 



+ (6ci 1 1 ) x + 34 2 1 ) x)(4-4 q )- (11) 



a 2 a = (6/tt) 



1/3 



,(2) 



, Alx -2fcTcxp(-44 2 lx /fcT) 
.(2) 



fcTexp(-6wXi x /fcT) , (13) 



where a is the fee lattice parameter. 

LTE thus allow to calculate all CNT input parameters 
from the knowledge of the solubility limit. The nucle- 
ation rate is then obtained from the equation 



J st (4) 



-167V. 



ag" uc (4) ; d x 

VkTa 2 a a 2 ' 



exp -- 



(a 2 a) 



3 fcT[AG nuc (4)] 2 



(14) 



4 




0.5 1 1.5 2 

%Zr 




0.5 

%Sc 



where Dx is the X impurity diffusion coefficient in Al. We 
thus obtain a fully analytical expression of the nucleation 
rate. Using the same experimental data, i.e. solubility 
limits and diffusion coefficients, as the ones used to fit 
the atomic diffusion model of kinetic Monte Carlo simu- 
lations, we can compare CNT predictions with nucleation 
rate observed in simulations^. A good agreement is ob- 
tained both for Al-Zr and Al-Sc binary alloys (Fig. [2]). 



The combination of LTE with CNT thus allows to build 
a quantitative modeling of nucleation relying on a very 
limited number of material parameters. Such a model 
can be directly applied to aluminum alloys where a LI2 
compound precipitates from a supersaturated solid so- 
lution as this is the case with Zr, Sc or other rare earth 
elements like Er, Tm, Yb and LiJ£. Li too precipitates in 
aluminum with a LI2 structure, but this system requires 
another statistical approximation than LTE. Indeed, this 
approach based on LTE, requires that the precipitating 
phase only slightly deviates from its perfect stoichiome- 
try and that the solute solubility remains low. Provided 
these conditions arc fulfilled, it could be applied to al- 
loys other than aluminum alloys. More generally, LTE 
demonstrate that the oscillating nature of the interac- 
tions in an alloy with an ordering tendency has to be 
taken into account by the CNT and requires a better 
statistical description than the Bragg- Williams approxi- 
mation which treats all interactions on the same footing. 



FIG. 2: Variation with nominal concentration and tempera- 
ture of the steady-state nucleation rate J st for AlaZr (top) 
and AI3SC (bottom) precipitations. Symbols correspond to 
kinetic Monte Carlo simulations and lines to classical nucle- 
ation theory 
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